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Abstract 

We describe exact kink soliton solutions to nonlinear partial differential equations in the 
generic form u t + P{u)u x + vu xx + Su xxx = A(u), with polynomial functions P(u) and 
A(u) of u — u(x,t), whose generality allows the identification with a number of relevant 
equations in physics. We emphasize the study of chirality of the solutions, and its relation 
with diffusion, dispersion, and nonlinear effects, as well as its dependence on the parity 
of the polynomials P{u) and A(u) with respect to the discrete symmetry u — > —u. We 
analyze two types of kink soliton solutions, which are also solutions to 1 + 1 dimensional 
6 4 and </> 6 field theories. 
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I Introduction 



Nonlinear differential equations are known Jj]] to describe a wide variety of phenomena not 
only in physics, in which applications extend over magnetofluid dynamics, water surface 
gravity waves, electromagnetic radiation reactions, and ion acoustic waves in plasmas, 
just to name a few, but also in biology and chemistry, among several other fields. In spite 
of the increasing development of mathematical techniques and concepts to solve nonlinear 
equations, exact solutions seem not to be the rule and numerical methods have been the 
most common approach to study their properties 0. 

Since the seminal work of Korteweg and de Vries (KdV) ||, in which a third-order 
nonlinear equation was studied to explain the shallow-water solitary wave experiments 
by Russel 0, solitary- wave or soliton solutions have been found in a number of nonlinear 
differential equations |T[, with phenomenologies related to liquid crystals, dynamics of 
growing interfaces and domain walls, eletromagnetism, nonlinear optics, acoustics, and 
elasticity, among others. These localized non- dispersive travelling waves can be of several 
distinct types 0, such as kink, pulse, breather, envelope, and dark solitons, all of them 
presenting the property that their shapes and velocities are preserved along propagation, 
and even upon collision with other solitary waves. In fact, they result from a precise 
balance among the competing elements that define the equation, namely the tendency 
of spreading due to the presence of a dispersive term and the action of nonlinear terms, 
which in general favour large amplitude disturbances and velocities that assure the sta- 
bility of the travelling waves with respect to small distortions in form. In some special 
cases, solitary waves do not support arbitrary speeds and can propagate only with defi- 
nite velocities determined from the parameters of the nonlinear equation. These so-called 
chiral solitons have received much attention lately due to their recent associations with 
the nonlinear Schrodinger equation and the fractional quantum Hall effect || |7|, || . More- 
over, we have also recently pointed |J to the presence of chiral kink soliton solutions in 
generalized KdV-Burgers-Huxley (gKdVBH) equations in the form 

u t + f x + 9xx - $ u xxx = h(u) , (1) 

where /, g, and h are smooth functions in u, and ut (u x ) stands for the partial derivative of 
u(x, t) with respect to time (position). The above equation combines dispersion, controlled 
by the real parameter 5, nonlinearity, described by the functions f(u), g(u), and h(u), 
and diffusion related to the second-order space derivative. 

This work is inspired in our previous study ||, although here we follow a different 
approach. We consider exact kink solutions to nonlinear differential equations of the 
generic form 

Ut + P{u) u x + v u xx + 5 u xxx = A(u) , (2) 
with polynomial functions defined as 

( 3 ) 

i=0 
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and 

N a 

A(u) = Y, a iU i - (4) 

8=0 

The general form of Eq. (2) allows the identification of several interesting cases 0. For 
instance, by considering propagating waves u(x,t) = u(x — ct) = u(y), with velocity c, 
it is easy to see that the gKdVBH equation is recovered from Eq. (2) for df fdu = P(u), 
g(u) = go + vu, 5 = —5, and h(u) = A{u). Furthermore, the standard KdVB equation 
corresponds to identifying P{u) = Xu, 5 = —5, and A(u) = 0, and the modified KdVB 
equation requires P(u) = 3Xu 2 , 5 = —5, and A(u) = 0, with the particular case v — 
accounting respectively for the standard and modified KdV equations. On the other hand, 
the Burgers-Huxley equation represents the situation in which P(u) = Xu, 5 = 0, and 
A(u) = h(u), with the case h(u) = corresponding to the standard Burgers equation, 
whose applications include nonlinear heat diffusion, shock waves, viscous effects in gas 
dynamics, and an important connection with the deterministic KPZ equation QlO | in one 
spatial dimension, known to provide the evolution of the profile of a growing interface or 
a domain wall of general nature. 

Other examples include the generalized Boussinesq equation jlL |5], 

u u ~ clu xx - p{u 2 ) xx - q{u 3 ) xx - hu xxxx = , (5) 

which by a trivial integration in y for travelling solutions u(x,t) = u(x — ct) = u(y) is 
related to Eq. (2) through the identifications P(u) = —c\ — 2pu — 3qu 2 , v = 0, 5 = —h, 
A(u) = a Q , where Co is the speed of the sound. The standard and modified Boussinesq 
equations correspond respectively to the cases q = and p = 0. They were first described 
in the context of the theory of long water waves, but also present innumerous applications 
in electromagnetism, plasma physics, elasticity, and other fields. Another generic equation 
commonly found in physics is [|l] 

2 dV 

u tt - a u xx + — = , (6) 
du 

which can also be compared to Eq. (2) by considering travelling waves with y = x — ct 
and identifying P(u) = c, v = c 2 — a 2 , 5 = 0, and A(u) = —(dV/du). As examples 
of Eq. (6), we recall the Klein-Gordon equation of field theory in which V(u) = fl 2 u 2 , 
also applied to describe standard vibrations for a displacement u with presence of an 
additional restoring force proportional to u, and the sine-Gordon equation corresponding 
to V(u) = (3 2 {1 — cosu), which embodies a wide range of applications including Josephson 
junctions in superconductors, dislocations in crystals, waves in ferromagnetic materials, 
laser pulses in two state media, and geometry of surfaces, just to name a few. Also, the 
standard nonlinear Schrodinger equation, 

%u t + uu xx = —u, (7) 
dp 
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corresponds to formally assigning for the travelling wave u(x — ct) the values P(u) = c—ic, 
v = p, 8 = 0, and A(u) = (dV/dp)u, where V = V{u*u) = V(p) is commonly expressed 
as a quadratic or cubic potential in the charge density p. 

Regarding the polynomial A(u) in Eq. (2), it first appeared as a generalization || 
of the study of monoatomic chains of equal masses interacting via a number of nearest- 
neighbor realistic potentials such as the Toda, Morse, and Lennard-Jonnes potentials flTT| . 
In the continuum spatial limit, the longitudinal displacement of the particle n from its 
equilibrium position, u n (t), is replaced by an analytic function u(x, t), and the series 
expansion of the interaction terms gives rise to the spatial derivatives of high orders. In 
this context, A(u) thus represents an on-site local interaction potential in the differential 
equation of motion that governs the system ||. 



In the following we generalize the procedure applied in Ref. [JT2J to some nonlinear 
equations, in order to investigate Eq. (2). The extension to equations with higher-order 
derivative terms, such as the fifth-order KdV-type equations [[|, is straightforward and 
will not be considered here. Our emphasis will be on the study of chirality of exact kink 
solutions to the general equation, and analysis of its relation with combined diffusion, 
dispersion, and nonlinear effects, as well as its dependence on the parity of the polynomials 
P(u) and A(u) with respect to the discrete symmetry u — > —u. In Sec. II we analyze 
the case in which kink solutions are of the type u(y) = atanh(Xay), whereas in Sec. Ill 
solutions in the form u(y) = {(a 2 /2)[l+tanh(Aa 2 ?/)]} 1//2 are considered. Finally, comments 
and conclusions are presented in Sec. IV. 



II Exact Kink Solutions of Type 



We start by considering travelling solutions u(x,t) = u(x — ct) = u(y) to the general 
third-order nonlinear differential equation in the presence of a local polynomial potential 
A(u), Eq. (2), 

r ,,du d 2 u ^d?u ., . , . 

[ - C + P(u)] Ty +V W +S W= A{u) - (8) 

From Eqs. (3) and (4), the orders of the polynomials, N p and iV a , are not arbitrary, 
but are instead strongly related to the characteristics of the solutions. For instance, if 
N p — N a — 2 one can search for exact kink soliton solutions in the form, 



since the property, 

leads to the general result, 



u(y) = atanh(Aay) , (9) 
^ = A(a 2 -« 2 ), (10) 



d n ii 

= <"> 
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Indeed, from Eq. (11) one sees that if N a > N p + 2 no solutions of the form presented in 



Eq. (9) can be found to Eq. (8). We have shown [|, |13[ that Eq. (9) also corresponds to 
kink solutions <fi = <p(y) in the context of the relativistic 4 field theory in 1 + 1 dimensions, 
which presents potential 

V(0) = iA 2 (0 2 -a 2 ) 2 , (12) 

that develops spontaneous symmetry breaking of the discrete Z 2 symmetry — > — 0. Here 
we mention that the kink solutions resembles the two asymmetric and degenerate vacua 
of this 4 system. 

The coefficients {p{\ and {a^} in Eqs. (3) and (4) are not arbitrary, but instead must 
obey a series of relations in order to allow Eq. (9) to be solution of Eq. (8). In fact, this 
is actually expected to occur to solitary-wave solutions, which can only exist for a precise 
balance of the diffusion, dispersion, and nonlinear effects. One can see that by explicitly 
substituting Eqs. (3), (4), and (9) in Eq. (8), to obtain 



Xa? 
X^ 

° 2 " P 2 + 85A 2 + ^1, (15) 



-c + p - 25X 2 a 2 , (13) 
Pi - 2vX , (14) 



Xa 2 " a? 

A 

T 



-Pi + 2vX , (16) 
-p 2 - 65 X 2 . (17) 



Although this technique can be easily generalized to any polynomial orders N p and N a , 
provided that in the case of solution given by Eq. (9) they obey the relation N p = N a — 2, 
we have restricted ourselves in Eqs. (13)-(17) to situations of physical interest, represent- 
ing quadratic nonlinearity in the first-order derivative term, N p — 2, as pointed from the 
examples listed in Sec. I. Indeed, considering the general N p > 2, N a > 4 case only gener- 
ates relations among highest-order coefficients {p{\ and {a{\ not involving the parameters 
of interest v, 5, and c. For instance, in this case Eqs. (16) and (17) would respectively 
read 



-pi + a 2 p 3 + 2z/A 



03 
A 

-1 = - P2 + a 2 p 4 - Q5X 2 , (19) 
A 

and new relations involving higher-order coefficients would emerge: 

y = -pi_ 2 + a 2 Pi , n > 4 . (20) 

Notice that Eq. (20) represents matches exclusively between coefficients of the polynomials 
P{u) and A{u) necessary in order to allow Eq. (9) as solutions to Eq. (8). No extra physical 
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information concerning diffusion or dispersion effects, or even chirality (see below) is 
gained from Eq. (20), contrarily to what happens to Eqs. (13)-(17), as we discuss below. 
Moreover, Eqs. (13)-(20) also indicate that the coefficients are not totally independent, 
but instead must obey a set of relations, which in the particular case of solutions given 
by Eq. (9) are given by, in the case of N a even, 

Na/2-l 

a 2 *a 2t+1 = 0, (21) 

and 

N a /2 

E fl2 V = . (22) 

i=0 

Consequently, setting the values of the N a /2 (N a /2 — 1) even (odd) Oj's automatically 
determines the values of the remaining even (odd) coefficients. 

We notice from Eqs. (13)-(20) that the diffusion and dispersion parameters, respec- 
tively v and 5, are related to definite and distict parities of the coefficients. For instance, 
v is associated only with odd coefficients, p±, p 3 , ai, 03, whereas 5 is related to the even 
ones, p , p 2 , Pa, a , a 2 , a 4 . No higher-order coefficients depend on v or 5. Indeed, this 
is consequence of properties of the solution considered, Eqs. (9) and (10), namely that 
each derivative multiplies the original function by a factor of order u, thus changing its 
parity with respect to the symmetry u — > —u, as one can also see from Eq. (11). Another 
important feature is that the velocity c of the kink solitons also depends only on even 
coefficients up to quadratic order, i.e., p , p 2 , a , a 2 . The fact that c is not arbitrary, be- 
ing instead determined by the parameters that define the nonlinear differential equation, 
indicates that solutions given in terms of Eq. (9) are actually chiral solitons. For instance, 
Eq. (13) can be used to determine the velocity of the localized travelling solutions, 

c = p -^(a Q + 25\ 3 a 4 ) . (23) 

We observe from the above equation that p plays the role of the velocity of the reference 
frame, as expected since it represents the w-independent term in the polynomial P(u). 
Furthermore, Eq. (23) also tells that the diffusion term associated with the parameter 
v does not play any role in fixing the velocity of the solitons, thus not contributing to 
their chirality. This is a consequence of the distinct parities of the coefficients associated 
with c and u, as commented above. Indeed, by considering the simplest case in which 
Po = do = 0, we point to the relevance of the presence of a dispersive media (5 7^ 0) in 
order to support chiral kink solitons of the type described by Eq. (9). Nevertheless, this 
does not seem to be a property shared by all types of kink solitons, as we comment on 
Sec. III. 

From the above discussions we also notice the remarkable fact that distinct choices of 
the parameters can lead to solitary- wave solutions with the same shape, Eq. (9), even if 
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they do not necessarily have the same velocity Indeed, this feature was first verified JTT 
with the finding that Toda, Morse, and Lennard- Jones solitons arising from distinct non- 
linear differential equations are very nearly alike, a result then attributed to the similarity 
in the shapes of their respective potential walls. 

Let us now consider some specific examples. The connection with the gKdVBH equa- 
tion, Eq. (1) with the ansatz u(x,t) = u(x — ct) = u(y), is easily established from the 
following identifications (see also Introduction): 

PC, = m , (24) 

vu = g(u) - g , (25) 
8 = -6 , (26) 
A{u) = h(u) , (27) 

where go is some unimportant constant. To illustrate with a previously reported || 
example of soliton solutions related to Eq. (1) and with form given by Eq. (9), let us 
consider the case in which 

P{u) = po+p 2 u 2 , (28) 
v = , (29) 
A(u) = ao + (I2U 2 , (30) 

that corresponds to the gKdVH equation with functions 

f(u) = fo+PoU + ^f, (31) 

9(u) = -g , (32) 
h(u) = a + a 2 u 2 . (33) 

Solving Eqs. (13)-(17) along with the functions expressed by Eqs. (28)-(30) leads to chiral 
solitons, Eq. (9), with velocity 



c = Po~ T. h sgn(A)a 2 \ . (34) 

3a 2 V V2 

The result reported in Ref. is restored by considering p = —2a 2 /3 2 , p 2 = 6j3 2 , a = a 2 /?, 
and a 2 = —(3. In this case the velocity of the chiral solitons are simply given by (5 < 0) 

c = sgn(A) v /Z 5". (35) 
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Ill Exact Kink Solutions of Type 

We now search for exact kink soliton solutions to Eq. (2) in the form 

u(y) = {(a 2 /2)[l + tanh(Aa 2 y)]} 1 / 2 . (36) 



As we have reported in Refs. || [T3|], differently from the solutions investigated in Sec. II, 
Eq. (35) corresponds to kink solitons of the </> 6 field theoretical system in the 1 + 1 
dimensional spacetime. In this case the potential is 

V(<f>) = ^A 2 2 (0 2 - a 2 ) 2 , (37) 

which also develops spontaneous symmetry breaking of the discrete Z 2 symmetry <fi — > —<p. 
Here, however, the kink solutions resembles both the symmetric and asymmetric phases 
this 4> 6 system engenders. 
In the present case we get 

^ = \u{a 2 - u 2 ) , (38) 
ay 

so that, by contrasting with the solutions studied in Sec. II, 

= u** 1 , (39) 

thus implying that each derivative operated preserves parity with respect to u — > — u by 
multiplying the original function by an even factor on u. As a consequence, one must 
now have N p = N a — 3 to allow Eq. (36) to be solutions to Eq. (8). In this case, however, 
since the third-order derivative of Eq. (36) gives rise to a polynomial of order u 7 , then 
the matching procedure involving physical ingredients such as diffusion and/or dispersion 
requires the presence of higher-order nonlinearity in A(u), N a = 7. By substituting 
Eq. (36) in Eq. (8), we obtain the following relations 

a = , (40) 



a 1 
Aa 2 

.2 



c + p + \a 2 v + \ 2 a^5 , (41) 



a 2 = \ a 2 pi , (42) 
c - p + a 2 p 2 - A\a 2 u - 13A 2 a 4 5 , (43) 

-Pi + a 2 p 3 , (44) 



A 

X 

A 

a 6 = -Ap 3 , (46) 

^ = - P4 - 15A 2 <5 . (47) 
A 



-p 2 + a 2 p 4 + 3Az/ + 27AV5 , (45) 
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First we notice that due to Eq. (39) all terms in the left-hand side of Eq. (8) depend on u, 
thus leading to the result of Eq. (40). Furthermore, as another consequence of Eq. (39), 
the diffusion and dispersion parameters, respectively v and 5, are both related only to even 
coefficients {pi} and odd coefficients {a^}, differently from the case considered in Sec. II 
in which they were associated with coefficients of definite but distinct parities (even for 
5, odd for v). Consequently, the physics related to diffusion and dispersion effects in 
Eq. (8) with solutions provided by Eq. (36) is irrespective to the presence of odd parity 
terms in the polynomial P{u), and even terms in A(u). Indeed, in this case the nonlinear 
differential equation, Eq. (8), presents odd symmetry with respect to u — > —u. The 
difference between the odd-symmetric and non-symmetric cases (even symmetry is not 
allowed due to the presence of the constant c) appears from the way nonlinearity enters 
the equation through P(u) and A{u). It is also interesting to observe that the velocity of 
the solitons c only depends on the even coefficients po and P2, and odd coefficients a± and 
az. This implies that, contrarily to the solutions discussed in Sec. II, both diffusion and 
dispersion effects are actually relevant to determine c: 

c = p - — (aj - i/A 2 a 4 - 5X 3 a e ) . (48) 
Act 

We also recall that, as for the previous case, the coefficients are not totally independent 
and relations given by Eqs. (21) and (22) still hold for soliton solutions in the form of 
Eq. (36). 

Finally, as an illustrative example we consider the case in which 

P{u) = po , (49) 

A(u) = -9AV<5w 3 + 24AV<5w 5 - 15A V . (50) 



and 
If we set 



v = -\a 2 5 , (51) 

then Eq. (8) with solutions given by Eq. (36) is equivalent to a gKdVBH equation, Eq. (1), 
defined by the functions 

/(«) = fo+PoU, (52) 

g(u) = g + uu , (53) 

h{u) = A{u) , (54) 

and with 5 = —5, if the velocity of the chiral kink soliton solution, Eq. (36), is 

c = Po , (55) 

which is also in agreement with Eq. (48) and the results reported in Ref. [[J. 
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IV Conclusion 



In this work we have investigated the presence of travelling kink soliton solutions u(x, i) = 
u{x — ci) to nonlinear partial differential equations of the generic form 

u t + P{u)u x + uu xx + 5u xxx = A(u) , 

with polynomial functions defined as 

N P N a 

These equations combine diffusion, dispersion, and nonlinearity in distinct ways, and 
include a number of examples of relevant differential equations in physics and other fields 
of nonlinear science. 

We have analyzed two types of kink soliton solutions, namely u(y) = atanh(Aay), 
which is also related to solutions of the relativistic 1 + 1 dimensional 4 field-theoretical 
system, and u(y) = {(a 2 /2)[l + tanh(Aa 2 |/)]} 1 / 2 , which is associated with kinks of the 6 
system. These kinks present distinct properties. For instance, while in the 4 system they 
connect the two asymmetric vacua, in the 6 model they relate the symmetric vacuum 
= to the asymmetric ones. 

The emphasis of our study was on the chirality of the localized soliton solutions, and 
its relation with diffusion, dispersion, and nonlinear effects, as well as its dependence 
on the parity of the polynomials P(u) and A(u) with respect to the discrete symmetry 
u — > —u. 
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